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A study is made of liow active membrane proteins can modify the long wavelength mechanics of 
fluid membranes. The activity of the proteins is modelled as disturbing the protein surroundings 
through non-local force distributions of which a force-dipole distribution is the simplest example. An 
analytic expression describing how the activity modifies the force-balance equation for the membrane 
surface is obtained in the form of a moment expansion of the force distribution. This expression 
allows for further studies of the consequences of the activity for non-planar membranes. In particular 
the active contributions to mechanical properties such as tension and bending moments become 
apparent. It is also explained how the activity can induce a hydrodynamic attraction between the 
active proteins in the membrane. 
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I. INTRODUCTION 

Biological membranes actively participate in many bi- 
ological processes. They constantly exchange material 
with their surroundings and they contain molecular ma- 
chines that consume free energy to perform different tasks 
0. Biological membranes are therefore not just simple 
self-assembled structures of lipids and proteins in thermal 
equilibrium, but should be considered as non-equilibrium 
systems. 

Micropipette experiments have been performed on 
lipid model membranes with membrane proteins that ac- 
tively pump ions across the membrane |3, 1^ ■ They have 
demonstrated that the activity of the pumps influence the 
mechanical properties of the membrane. The membrane 
protein investigated in 0, bacteriorhodopsin, was a light- 
driven proton pump, while the protein studied in .Sj was 
the ATP-driven calcium pump Ca^+-ATPase. In both 
micropipette experiments it was found that the mem- 
brane tension was less sensitive to changes in the visible 
area of the membrane when the pumps were active. This 
was interpreted as the shape fluctuations of the mem- 
branes being enhanced by the activity. Therefore the re- 
sults were presented as increased effective temperatures, 
which could be about two to three times the values found 
when the micropipette experiments were performed with- 
out the pumps being activated. 

In 3' theoretical model was proposed which explains 
the enhancement of the fluctuations as being a result of 
the active proteins pushing on their surroundings. The 
pushing was modelled in the simplest possible way, as a 
force-dipole acting on the fluids surrounding the mem- 
brane. A force monopole was excluded because it would 
require an external source of force in the system. Com- 
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pleting the model by setting up stochastic equations of 
motion for the membrane shape to first order in devi- 
ations from a planar membrane, including an averaged 
contribution of the active force-dipoles together with 
thermal noise in the description, it was shown that this 
model provides a possible explanation of the enhance- 
ment of the fluctuations. However, as noted in the 
non-thermal noise from the activity may also contribute 
to the enhancement of the fluctuations. 

To test the force-dipole model more stringently it 
would be necessary to perform other experiments on 
these active membrane systems. One possibility would 
be to measure the fluctuation spectrum directly using 
video microscopy as proposed in Another possibility 
would be to study how the activity modifies the overall 
average shape of the membrane. For lipid membranes in 
equilibrium this shape is well studied both theoretically 
and experimentally, see for instance for a review. A 
way to make the average shape sensitive to switching 
on activity would be to prepare the membrane in a state 
where a small change in control parameters, such as tem- 
perature or osmolarity, will result in a sharp transition 
in shape, for instance from a prolate to an oblate (24l |. 

To find the predictions of the force-dipole model for 
such experiments it would be beneficial to be able to 
work analytically with the model for membrane confor- 
mations that are not close to being planar. However, as 
the force-dipole model was formulated in 4], this is not 
easy to do directly, because the equations of motion for 
the surrounding bulk fluids would then have to be solved 
in these non-planar geometries including the presence of 
the force-dipoles. 

The purpose of the present paper is to show how this 
difficulty can be avoided by using a more indirect ap- 
proach where the formalism developed in is applied. 
The idea behind this approach is to develop an ideal- 
ized mathematical formulation of the force-dipole model 
where all effects related to the presence of the membrane 
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are assigned to an infinitely thin surface. The space 
around this surface is filled with the bulk fluids and ex- 
cess quantities are then attributed to the surface such 
that the amounts of extensive quantities like energy, mo- 
mentum, number of molecules etc. are conserved. This 
idealized formulation will also be called the Gibbs formu- 
lation of the model, since the philosophy behind it, ap- 
plied to a non-equilibrium situation here, is the one that 
Gibbs developed to treat the thermodynamics of inter- 
faces |9|. For the force-dipole model the most important 
adjustment needed to arrive at its Gibbs formulation is 
to replace the non-local contribution of the force-dipoles 
in the bulk fluid momentum conservation laws with an 
excess current of momentum on the membrane surface. 

The outline of the presentation of this approach is as 
follows. We will give a brief review of differential geom- 
etry of surfaces in section ^1 to establish the notation. 
Then we will continue with formulating a generalized 
version of the force-dipole model of ^] in section IIIII 
This formulation of the model will be called the semi- 
microscopic formulation. Microscopic, because molecu- 
lar length scales enter directly, namely the lengths over 
which the forces are distributed. Semi is added in front 
because, apart from these activity induced forces, the 
membrane is treated as being infinitely thin. In section 
II VI we then show how an expansion can be performed in 
these molecular length scales divided by the curvature 
radii of the membrane to arrive at the Gibbs formula- 
tion for the model. In the Gibbs formulation some of 
the physical consequences of the force-dipole model be- 
come clearer than in the semi-microscopic formulation. 
One of these consequences is an influence of the activ- 
ity on the tension of the membrane, which is discussed 
in section Another consequence of the force-dipole 
model, discussed in section IVTl is a hydrodynamic inter- 
action force between the membrane proteins induced by 
the activity, offering a possible explanation of two exper- 
iments reported in 10] on the clustering and diffusion 
of active bacteriorhodopsin molecules. A conclusion is 
given in section IVIII and finally an appendix is added 
where the bulk hydrodynamics is solved directly for the 
semi-microscopic formulation in the case of a nearly pla- 
nar membrane. Comparing this solution with the Gibbs 
formulation derived in section HVI it can be seen directly 
that the Gibbs formulation and semi-microscopic formu- 
lation are equivalent in this case. 

The Gibbs formulation derived here will be applied in 
another paper [TH . where the fluctuation spectrum of 
a quasi-spherical vesicle with active proteins acting as 
force-dipoles is derived. 



II. BRIEF REVIEW OF DIFFERENTIAL 
GEOMETRY OF SURFACES 

In this section we briefly review the mathematical lan- 
guage of two-dimensional differential geometry, which 
will be used throughout the rest of the paper. 



The dynamic shape of the surface is represented by a 
space- vector function R = where the vari- 

ables f ^ and f ^ parametrize the surface and t represents 
time. At each point on the membrane surface we have 
a basis for three-dimensional vectors consisting of two 
tangential vectors 



la - Uatt - , 



(1) 



where a = 1, 2, and a unit vector normal to the surface, 

1*1 Xt2| ^ ^ 

The metric tensor of the surface is defined by 

5q/3 =ta-ti3. (3) 

It has an inverse, g"^ , which satisfies 

g'^f'g^^ = <5^ , (4) 

where i5" is the Kronecker delta and where the repeated 
Greek superscript-subscript indices imply summation fol- 
lowing the Einstein summation convention. The metric 
tensor and its inverse are used to raise and lower Greek 
indices as in the following example: 



t = g ^ti3 , ta = 5a/3C 



The curvature tensor Kaf3 is 

Ka0 = n ■ dadpR 



(5) 



(6) 



From Kap the scalar mean curvature H and Gaussian 
curvature K can be obtained: 
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H = -g'^'^Kap , 
K = det g°''^Kp^ 



(7) 
(8) 



The expression for covariant differentiation of for in- 
stance a vector w = w"ta is given by 

Dawl" = dawf" + wT^^^ , (9) 

where the Christoffel symbols are defined as 

= ^9'^^ {dpgSa + da9f35 " dsgafi) ■ (10) 

Finally, the area of a local differential element of the 
surface is given by 

dA = ^dede , (11) 

where g = det gaf} is the determinant of the metric tensor. 
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III. SEMI-MICROSCOPIC FORMULATION 

The force-dipole model proposed in 3 states that the 
important contribution of the activity of the active mem- 
brane proteins to the mechanics and dynamics of the 
membrane shape is that the proteins push on the sur- 
rounding bulk fluids. Mathematically, this was formu- 
lated by the addition of a source of force .Fact in the 
equations of motion of the bulk fluids. Due to the mi- 
croscopic size of lipid-protcin membranes we will assume 
that we are at low Reynolds number where these equa- 
tions become 



r,V^v± - Vp± + Ft, = , 



supplemented by the incompressibility condition 
V • -Ui = . 



(12) 



(13) 



Here v± = v±{r, t) are the velocities of the bulk fluid at 
position r and time t on the side that n points into (+) 
and the other side (— ) of the membrane, p± — p±{r,t) 
are the corresponding pressures, -F^ct = ^acti'''^ ^) ^^'^ the 
appropriate restrictions of the active force density .Fact 
to the ±-side of the membrane and 77 is the viscosity of 
the fluids. 

The speciflc expression for i^act will be generalized 
slightly here compared with jj]. We will write it as 



Fact(r) 



M 



dA / dh F^ctnS-" {r- {R+ hn)) . (14) 



The first integral in Equation (|14|l is over the area of 
the membrane M, and the second integral is along the 
normal direction to the membrane with h being the dis- 
tance to the membrane surface. Together with the delta 
function these integrals can be thought of as a conver- 
sion of coordinates from the Cartesian position vector r 
to the membrane related coordinates (C^, C^, h). F^a will 
be taken to be an arbitrary function of h and fields on the 
membrane, say -Fact = Fi^ctih, H,np+ ,np-), and it there- 
fore depends implicitly on (C^,^^) and t. np± represent 
area densities of the active proteins in the membrane, 
with -I- and — indicating the two possible orientations of 
an asymmetric transmembrane protein which is incorpo- 
rated in a membrane. 

We can at any point specialize to the more specific 
force-dipole model of 0|| by taking 

fact = (FaUA + 2i/F>s) [S {h - w^) ~5{h + w^)] , 

(15) 

where Fa and F'^ are constants representing the strength 
of the active forces and their curvature dependence, 
and are constant lengths giving the distances from the 
membrane where the forces act. For the protein densities 
we have introduced ns = np+ -t-Up- and ua = 'T-p+ —np- . 

At this point it should be pointed out that there is a 
problem with Equation H15|l : it is not invariant with re- 
spect to exchanging the two sides of the membrane (see 



«,,+ protein 



i 









Eq. (15) 



protein 



Eq. (16) 



FIG. 1: Illustration of the points of action and directions of 
the forces in the force distributions of Equations H15|l and 



Figure If locally the physical situation is symmetric 
with respect to this exchange, i.e. tia = and the mem- 
brane is planar, then the equation should be invariant 
with respect to exchanging the two sides. As also noted 
in 12], this problem can be resolved by using instead the 
symmetric dipole model 



{Fa + 2HF'^np+ [5{h-w^] 

+ {Fa - 2HF'a) Up- [5{h~v 



- 5[h + w^)] 
^) ~S{h + w'^] 



(16) 



In section we will discuss what the difference is in the 
consequences of the two specific models. Equations (|15|) 
and H16|l. for the mechanics of the membrane. 

In the Navier-Stokes equation 1)12(1 i^act acts as a local 
source term for momentum. Globally, however, there can 
be no production of momentum. Or equivalently, formu- 
lated as in Newtons third law, for every force there has 
to be an equal but opposite counterforce. This condition 
will be ensured by assuming that 



dh Fact = 



(17) 



To complete the description of the bulk hydrodynam- 
ics we need to supplement Equations H12() and (|13|l with 
boundary conditions on the velocities v± at the mem- 
brane position R and at the boundaries of the full 
membrane-bulk fluid system. We will not choose any 
specific boundary conditions here, but for later use we 
will introduce the following notation for the boundary 
conditions at the membrane 



v±\ 



-R 



(18) 



where U± = U± are unspecified. 

Equations (|12|1 and (|13|l can be solved for the case 
of an almost planar membrane, expanding to first or- 
der in deviations of the membrane shape from a plane. 
This is carried out in the appendix of this paper. It 
was also carried out in a different manner in where 
the result was used to explain data obtained from micro- 
pipette experiments on active membranes. However, the 
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large scale geometry of the membranes in micropipette 
experiments are not planar. Treating them as planar is 
an approximation justified for studying membrane un- 
dulations with short wavelengths compared to the typi- 
cal length scales of the overall membrane geometry. For 
the case of micropipette experiments it is sufficient to 
know the behavior of the short wavelength undulations 
because the long wavelength undulations are suppressed 
by the tension that is induced when the membrane is 
aspirated to the pipette. But as explained in the Intro- 
duction there are situations where it is useful to know the 
behavior of the long wavelength undulations or the over- 
all geometry of the membrane shape. The problem with 
the semi-microscopic formulation is that directly solving 
the bulk fluid hydrodynamics analytically for this formu- 
lation seems to be quite difficult, if not impossible, for 
a typical membrane geometry like one being close to a 
sphere. In the next section we will show how this prob- 
lem can be circumvented by a more indirect approach to 
the problem. 



The free energy for a lipid membrane should include 
terms representing bending resistance and tension 



F ^ dA [2kH^ + (To + • • • ] 
Jm 



(22) 



where k is the bending rigidity, ctq a tension and the dots 
represent terms involving free energy variations with re- 
spect to changes in other fields such as the density fields 
np± . The force /^j^ is associated with internal dissipation 
in the membrane. However, this dissipative force is often 
taken to be negligible in comparison with the dissipation 
in the bulk fiuids because of small thickness of the mem- 
brane. Instead dissipation enters through T^, which are 
the forces from the bulk fluids on the membrane. Given 
that the stress tensors of the bulk fluids are 

^-p^\ + rj\Vv± + {Vv±f] , (23) 

where I is the identity tensor, we can flnd as 

= ±n • T± . (24) 



IV. GIBBS FORMULATION 

In this section we will derive the Gibbs formulation 
corresponding to the semi-microscopic formulation of the 
generalized force-dipole model of the last section. 

The designation Gibbs formulation will in this context 
be taken to mean a mathematical model of the membrane 
system in which the membrane is treated entirely as be- 
ing a surface with zero thickness. The semi-microscopic 
formulation of the generalized force-dipole model in the 
last section is not a Gibbs formulation, because the active 
membrane proteins are acting with forces in the bulk flu- 
ids across a finite distance transverse to the membrane. 
Thus in the corresponding Gibbs formulation these active 
forces will have to be absent from the bulk fluid equations 
of motion, i.e. instead of Equation (|12|l we would have 



r]V^v± - Vp^ = 



(19) 



and the effect of the active forces will have to be present 
in the equations of motion for the membrane surface and 
the boundary conditions for the bulk fluids instead. 

One of the equations of motion for the membrane sur- 
face is the conservation of excess momentum. At low 
Reynolds number where inertial and convective terms can 
be discarded this law becomes a force balance equation, 
which can be written 1251 



/act = . 



(20) 



Here is the elastic restoring force of the membrane 
derivable from the membrane excess free energy F by 
functional differentiation Mi 



fr 



1 SF 



(21) 



The main character in this paper is /^ct; which is the 
force on the membrane induced by the activity of the 
membrane proteins. Note that the force balance equation 
for the membrane in the semi-microscopic formulation 
would be Equation 1^0)1 without /act: because in this 
formulation the activity induced force is included through 
. Thus in the semi- microscopic formulation the forces 
acting directly on the membrane is the same as those 
included in earlier formalistic works like for instance Il3l 

Q. 

Figuring out what the force /^ct is m the Gibbs for- 
mulation is, however, not trivial. A priori /^ct can be 
divided into two contributions 



/act ^ ""a 



^Q-*- ac- 



(25) 



where cr act is the excess source of force per area due to 
the activity and T"^^ is the excess stress induced by the 
activity. <Tact can be found by projecting the volume in- 
tegral of i^act down to an area integral on the membrane 
surface giving 



dh (1 - 2hH + h^K) Fa 



(26) 



since the volume element can be written dV — 
(l — 2hH + h^K) dAdh. However, as a consequence of 
Equation (|17() . there is no excess source of momentum 







(27) 



This leaves us with T^di which can be found by requir- 
ing that the stress on surfaces intersecting the membrane 
orthogonally is the same in the semi-microscopic formu- 
lation and the Gibbs formulation. This implies 

T:,, = jdh [g'^P - h {2Hg'^^ - K'^^)] tp ■ Tact , (28) 
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Pressure Pg 

- ■ ^''"^ \^^n{l + 2hH) F^tdh 



-r"^'"'-; '\ Pressure Pq" — (1 + 2hH)E^f-tdh 




'■■ 'j^^^^ Area c=l A ^ ^ 


h 



FIG. 2: Illustration of considerations leading to Equation 
l|30|l . Note that H is negative. 

where Tact is the bulk fluid stress induced by the activ- 
ity of the membrane proteins. More precisely, Tact is the 
bulk fluid stress tensor in the semi-microscopic formula- 
tion minus the bulk fluid stress tensor in the correspond- 
ing Gibbs formulation, with the boundary conditions for 
Equations H12(l and H19(l chosen such that the solutions 
for the bulk fluid motion are identical (for identical mem- 
brane shapes and protein density fields) in the two for- 
mulations at distances from the membrane greater than 
those connected with the activity induced forces Q. 

It may seem that to find T^ct as a function of the fields 
on the membrane we still need to work out the hydrody- 
namics in the semi-microscopic formulation completely. 
However, if we do not insist on obtaining the exact for- 
mula for Ta^.^, we can obtain an answer as a moment 
expansion of -Fact in the distance h from the membrane 
in a simpler way. 

The strategy we will follow is to work out the formula 
for the moment expansion of T^ct up to the second mo- 
ment in a special case, and then use a combination of 
symmetry arguments and dimensional analysis to argue 
that the formula obtained in the special case is actually 
the general formula. We will also use the linearity of 
the problem to divide the force i^act into infinitesimal 
contributions of the form 

F,ct,h ^dh ( dA F,,tnS^ {r-{R+ hn)) , (29) 
Jm 

where dh should be understood as an infinitesimal length. 
Then we can find the corresponding contribution T"^^ ^ 
to Tact fo'" each of these, and the full result is obtained 
by performing an integral over h in the end. 

The special case we will study is that of a homoge- 
neous membrane with constant curvature, such that n^i , 
H and K are constant and D^K^p — 0. Examples of 
such geometrical surfaces would be planes, cylinders and 
spheres. The boundary condition for the Navicr-Stokes 
Equation l|12|l is chosen to be that the velocity should 
vanish at the position of the membrane. 

The solution to the hydrodynamics in this special case 
will simply be that the bulk fluid will not move, v± = 0, 
but there will be a jump in the pressure at a distance h 



from the membrane to balance the active force. The dis- 
continuous jump for a planar membrane will be F^ctdh, 
with the lower pressure close to the membrane if -FkctJ^ 
points away from the membrane. If the membrane is 
curved with constant curvature, then the force acting at 
the distance h will be changed slightly, due to the fact 
that the area of the surface at the distance h differs from 
the membrane area (see Figure^}. To first order in h we 
will have a jump in the pressure which is given by 

A ± _ / T{l + 2Hh + 0{h^)) Factd/i , ±h>0, 
^^'^ " \ ±h<0. 

(30) 

This is the change in pressure between the membrane 
and the distance h induced by the active force i^act./i- 
From Equation H28|l we have that the corresponding ex- 
cess stress, Tact./i is 

T:,t.H = / dh' (g"^ - h' (2i75"^ - K'^^)) tp ■ Tact,. , 

(31) 

where Tact,/i is equal to —^Pf^ times the identity for val- 
ues of h' between and h, and zero otherwise. Perform- 
ing the integration over h' we get to second order in h 
for the tangential part T^",f^ = T^.^ • t^, 

T^lk = F-cth dhg^f" + (Hg'-P + ^if"^^ F^^.h" dh . 

(32) 

Finally, performing the integral over h we arrive at 

Ta"? = ■ = adipff"'^ + [ng'^f' + \k"P^ Q , (33) 
where 

CTdip = j '^^ hFact , (34) 
Q = J dhh^F,,t , (35) 

are the dipole and quadrupole moment of the force dis- 
tribution Fact. 

Even though formula (|33|l was derived for homoge- 
neous membranes with constant curvature, it has to be 
true for any geometry and distribution of active protein 
molecules, as long as the geometry and distribution are 
smooth enough such that we can expand T"Jl in curva- 
ture and gradient operators. This is true because other 
contributions will have to involve the gradient operators. 
And since we are building a second-rank tensor, there 
will have to be at least two of them (the only vectors we 
can use are the gradients of the force distribution and the 
curvature). From dimensional analysis we can then see 
that the lowest order quantity we can build using gra- 
dient operators involves the gradient used twice on the 
third moment of Fact, which is one order higher than 
what we are expanding to. 

The above symmetry argument also tells us that we 
cannot derive the normal component T,^act = ^act ' 
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directly by using a homogeneous membrane, since there 
wiU have to be at least one gradient in the formula for 
T^n a.ct to end up with a surface vector. However, there 
is another way, which goes via obtaining a general for- 
mula for the bending moments, M"^, arising due to 
the activity. These bending moments are defined such 
R X T"^.^. + M^J^n X tp is minus the non-convective part 
of the excess angular momentum flow induced by the 
activity, i.e. M"f^ is related to the activity induced inter- 
nal angular momentum current in the Gibbs formulation. 
By requiring that the angular momentum flow through 
surfaces intersecting the membrane orthogonally is the 
same in the semi-microscopic formulation and the Gibbs 
formulation one can derive that 15] 

M:f, = J dh' h' is"'' ~ h' (2i/.g«^ - if "^)) • Taet • ■ 

Since there is no excess source of torque in the system, 
we will have that the internal torque in the membrane 
vanishes. The excess internal torque is the excess torque 
on the membrane minus the torque of the excess force on 
the membrane [l^ . i.e. its activity induced part is 

(r X TZ,, + AC^n X t^) - X f . (37) 

The balancing of torques in both the semi-microscopic 
formulation and the Gibbs formulation implies that the 
activity induced part of the internal torque vanishes by 
itself. Equating the tangential part of Equation l(S7|) by 
zero we find that the bending moments are related to the 
normal component of the stress by the formula 

T„",,t = DpM:!, . (38) 

Inserting the solution for Tact.h for the homogeneous 
membrane case in Equation H36|) and integrating with 
respect to both h and h' , we obtain 

m:1 = \Qa-^ . (39) 

Since the bending moments again form a second rank 
tensor we can use the same symmetry argument as for 
T^t to argue that is the general formula for Af^^.^, 
and we thus have the desired formula 

r,:act = \q^Q , (40) 

correct to the second moment of i^act for any smooth 
membrane composition and shape. 

As briefly touched upon after Equation H28() . we need 
to have a boundary condition for the bulk fluid in the 
Gibbs formulation that ensures that we get the same be- 
havior of the bulk fluid far away from the membrane as 
in the semi- microscopic formulation |^. 

To find the boundary condition in the Gibbs formu- 
lation, we will again study a special case in the semi- 
microscopic formulation. This time the special case will 



be a planar membrane positioned at ^ = in a Cartesian 
coordinate system, {x,y,z), but with an inhomogeneous 
distribution of active proteins, such that the force obeys 
Fg^tdh = Cx with C being independent of x and y. For 
the case of a boundary condition where the velocity van- 
ishes at the membrane we can now solve Equation H12|) 
when i^act is restricted to Fa,ct,h- For the upper part of 
the bulk fluid a solution is that p'^ and are constant 
outside h and for z < h they are 

= constant — Cx , (41) 
C 

v+ = z{z-2h)x. (42) 

2r) 

This implies, due to the continuity of at z = ft,, that 
in the rest of the bulk we will have 

v+ = ^h'x . (43) 

The parametrization invariant boundary condition in 
the Gibbs formulation that agrees with Equations H43|) 
and is 

v±\r=R^U±±^t"D^Q^ , (44) 

where 

Q± = ± / dh ft^Faet . (45) 

Jo 

Again we can use symmetry arguments and dimensional 
analysis to justify that Equation is the general for- 
mula for the boundary condition up to the second mo- 
ment of i^act true for any smooth configuration of the 
membrane. 

Even though the modification of the boundary con- 
dition in Equation (|44|l is a strict consequence of the 
generalized force-dipole model formulated in section ITTll 
one can argue that for some proteins it is an anomaly 
arising because the model is too simplistic. In the case of 
for instance bacteriorhodopsin, we have that the physical 
extension of the proteins outside the width of the lipid 
membrane is limited ,16]. In this case it seems reasonable 
that an improvement of the Gibbs formulation would be 
to discard the modification in Equation (|44l) and use the 
original Equation 118|) instead. Whether or not the same 
argument can be used to question how realistic the mod- 
ification of the membrane force-balance coming from the 
force-dipole model, /act in Equation (|20|) , is relies on how 
close the mechanical properties of the lipids in the mem- 
brane resemble those of an incompressible fiuid. We will 
not try to resolve this here. 

It was mentioned at the end of section IIIII that bulk 
fluid hydrodynamics in the semi-microscopic formulation 
can be solved directly for an almost planar membrane. 
As a check of the correctness of the Gibbs formulation 
derived in this section, this solution is presented in ap- 
pendix^ The resulting modification of the force balance 
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equation for the membrane that the active proteins are 
found to induce are in complete agreement with the re- 
sults derived in this section. 

The momentum conservation law is not necessarily the 
only conservation law in the Gibbs formulation that is 
modified through the activity of the proteins. A pos- 
sible mechanism through which the active proteins can 
attract each other, leading to modification of their diffu- 
sion equations, is discussed in section IVTl 



V. MODIFICATION OF THE TENSION 

Some of the consequences of the force-dipole model for 
the mechanics of the membrane emerges readily in the 
Gibbs formulation. In this section we will discuss the 
effect of the dipole moment ddip on the tension of the 
membrane. 

From Equations H33|l and H4(J|) we see that tTdip only 
contributes to the isotropic part of the tangential stress 
in the membrane. This part is a two-dimensional analog 
of pressure in three dimensions 0,0, 'l8j, but due to 
the difference in sign it is called tension. Thus the role of 
the dipole moment (Tdip is exactly that of adding a local 
contribution to the tension of the membrane. 

For the specific dipole model of |3| , Equation H15() , the 
active contribution is 



CTdip 



- (.J 



(Fan A + 2i?F>s) 



(46) 



If TT-A = and the membrane is planar, which they ac- 
tually go to some length to make certain are reasonable 
assumptions in the micropipette experiments of , then 
the prediction of Equation (|4(j|l is that CTdip = on av- 



erage. However, as mentioned in section Hill the specific 
model of Equation H15|l is not realistic because it does 
not respect the symmetry of exchanging the two sides 
of the membrane, which should be there when riA = 
and H — 0. Therefore we will instead use the expression 
for (Tdip that one obtains from the symmetric model of 
Equation which is 

(47) 



CTdip = {w^ + w^) (FaU^ + 2ffF>A) 



It should be mentioned that even though the expres- 
sion for the active contribution to the tension is different 
for the two dipole models, their predictions for the ef- 
fective temperature in the micropipette experiments are 
the same. This is the case because the tension is fixed 
in the micropipette experiments by other control param- 
eters such as pressures, and because the expression for 
the other parameter contributing to the mechanics, the 
quadrupole moment: 



(Fan A + 2i?F>s) , (48) 



is identical for the two models. 

To give an estimate of the active contribution to the 
tension, Eauatio nl[47|l . we will use an estimate of a quan- 
tity "Pa made in which was found to agree well with 



the micropipette experiments on bacteriorhodopsin. The 
quantity is defined as 



'a — 



2w 



(49) 



where w ~ 5 nni is the membrane thickness, and the 
estimate is 

Va ^ lOkBT . (50) 

Making the assumption — ~ w/2 we get 

Fa (w^ + w^) = , ~ AOkBT ~ 1.6 X 10-^^.1 . 

w' — 

(51) 

In 01 they achieve concentrations around ns = 10^^ m^^. 
Multiplying this number with the one from Equation l|51|) 
we find an estimate for the tension 



CTdii 



Fa {w^ + w^) ns ~ 1.6 • 10"^N/m 



(52) 



If the membrane is in a floppy state with a large ex- 
cess of area stored in the fluctuations of its shape, which 
would for instance be the case if the fluctuations of the 
membrane were to be measured in a video microscopy 
experiment, then the membrane can contract its area in 
response to the increased tension. Thus the true increase 
in tension will not be as high as indicated by Equation 
(|52ll in that case. Assuming that almost all of the ac- 
tive contribution to the tension is canceled by a con- 
traction we can estimate the relative change in area by 
Aa — ~(Jdip I Ka, where Ka is the membrane expansion 
modulus. A typical value for Ka would be Ka = 0.2 N/m 
[T9l |. Using the tension from 1)5 2|l . we then get a relative 
decrease in area of around 0.8 percent. To estimate the 
the small part of the active contribution to the tension 
which is not canceled we assume that we have a nearly 
planar membrane sitting on a frame of area ^o- Expand- 
ing the free energy in Equation H22|) . assuming that the 
flelds contributing the terms indicated by dots have al- 
ready been integrated out, to second order in deviations 
e = e{x, y) from a planar shape one can find using the 
equipartition theorem that the fluctuation spectrum of 
the membrane is 

(KP)=^o (53) 



where = J dxdy e^^'^^^^'^yy^e{x,y) and q — ^Jo^ 



The excess of area stored in these fluctuations can then 
be calculated as 



8ttk 



ksT (To 
'7, log 



2 



(54) 



where V|| — t"da and gmax is a cut-off at high wavenum- 
bers where the free energy in Equation (|22|l is not ex- 
pected to be valid anymore. A cut-off at small wavenum- 
bers is not necessary due to the suppression of the fluc- 
tuations by the tension there. Equation H54|) can be used 
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to obtain k from micropipette experiments by measur- 
ing the relative excess area as the tension is changed. 
In |2| this was used to obtain the value k ~ lOkBT for 
the bending rigidity of a passive lipid-bacteriorhodopsin 
membrane. When the bacteriorhodopsin molecules were 
activated by shinning green- yellow light on the membrane 
then it was found that the change in relative excess area 
with tension now behaved as 



log 

UTTK CTinitial 



(55) 



where the effective temperature introduced here was 
measured to be T"^^ ~ 2T. Isolating the tensions we 
find from Equation (|55l) that a 0.8 percent decrease in 
the excess area corresponds to an increase in tension by 
a factor 



gfinal 
f-^initial 



exp (0.008 •87r-10/2)~3. 



(56) 



Thus if the initial tension is much smaller than dtjip then 
the assumption that almost all of the active contribution 
is canceled by the contraction of the membrane area is 
not violated by this factor three increase of the initial 
tension. 

Note that the assumption — u'-'' ~ u;/2 is a con- 
servative estimate, meaning that if the dipole model in 
Equation Hlt)|) is appropriate and if Va has the value given 
in Equation (|50|l . then the estimate in Equation H51|l is 
a lower bound on the dipole moment of the individual 
active proteins. However, if the lower bound is contra- 
dicted by experiment, the force-dipole in Equation 
could be replaced with a more complex force distribution 
giving a smaller dipole moment but the same quadrupole 
moment. 



VI. HYDRODYNAMIC INTERACTION 
BETWEEN ACTIVE PROTEINS 

As mentioned in the introduction, two experiments on 
the diffusion dynamics and the interactions of bacteri- 
orhodopsin molecules in their passive and active states 
were reported in p^ . The experiments showed that the 
diffusion of the bacteriorhodopsin molecules slows down 
when they are activated, and that the bacteriorhodopsin 
molecules tend to cluster more in their active state. Both 
of these effects can be explained by an attractive force 
between the bacteriorhodopsin molecules arising in the 
active state. 

A possible mechanism through which this attraction 
can take place is a coupling between bulk fluid pressure 
and diffusion currents in the membrane which was also 
discussed in One way to state this mechanism is 
that there should be a term in the chemical potential 
of the membrane proteins which is the pressure of the 
surrounding bulk fluid times the volume of the protein. 
Thus when the activity of a protein lowers the pressure in 
its surroundings. Equation (|30|l . the chemical potential of 



other molecules in its vicinity will be lowered and there- 
fore there will be an attractive force toward the active 
protein through this mechanism. However, the lipids in 
the membrane will also be attracted, and thus a counter- 
pressure will arise inside the membrane opposing the at- 
traction. But if the proteins stick out of the membrane 
they will still feel a net attractive potential equal to the 
lowering of the surrounding pressure times the volume of 
the part of the protein which is outside the average lipid 
thickness of the membrane. 

To give a rough estimate of the importance of this 
attractive mechanism, we can compare it with the en- 
tropic repulsion of the proteins. Ignoring that proteins 
can insert themselves with two possible orientations in 
the membrane we can take the chemical potential of the 
proteins to be 



^kBT\ognp+pAV , 



(57) 



where AV^ is the volume of the part of the protein which 
is outside the lipid part of the membrane and p is the 
pressure of the bulk fluid next to the membrane. Taking 
the membrane to be planar we get from Equations Ij^Ufl 
and (|16|l that p — po — FaUp where po is the pressure 
in the absence of activity, which we will take to be a 
constant. Working to first order in deviations from the 
homogeneous distribution of proteins, we can write the 
constitutive relation for diffusion as Q 

-g^^D^ (rip^VP) ~ DpD^d'^np , (58) 
where flp is a kinetic coefficient for protein diffusion and 



Dp = ftp 



-FaAV 



(59) 



is the corresponding diffusion constant. Thus to esti- 
mate the importance of the active contribution relative 
to the entropic, we can compare the two terms in the 
parenthesis in Equation H59fl . Reusing the protein den- 
sity Up ~ 10^^ m~^ we find 



~ 100 nm^ • kBT . 



(60) 



According to 'l(| the parts of a bacteriorhodopsin 
molecule extending into the aqueous surroundings are 
limited. Let us, however, say that these parts extend ap- 
proximately 1 nm out. Assuming that the cross-sectional 
area is approximately (5 nm)^ we get a volume of ap- 
proximately AV ~ 25nm'^. To estimate Fa we can use 
the estimate in Equation (|51|l and set w'^ + u 
Doing this we get 



5 nm. 



FaAV 



40 fesT 



• 25 nm'* = 200 nm^ ■ k^T . 



(61) 



5 nm 

The lesson that we can take from the above estimate 
is that this force-dipole induced hydrodynamic attrac- 
tion can matter for the clustering and diffusion of bac- 
teriorhodopsin molecules. Working out the importance 
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of this mechanism more precisely is difficult, because it 
would require knowledge of what additional forces, for 
instance membrane induced interactions , there is be- 
tween the proteins besides the active and entropic contri- 
butions. In the case of bacteriorhodopsin there has to be 
additional attractive forces, since it is known that bacte- 
riorhodopsin molecules also form clusters in their passive 
state 10, 21]. 

It should be mentioned that there is also mechanisms 
that can actively enhance diffusion. Rotating motors for 
instance enhance the fluctuations in the local flow around 
them creating an active contribution to diffusion j23 |. 



VII. CONCLUSION 

In this paper the force-dipole model of 3 was gener- 
alized to arbitrary distributions of forces in the direction 
transverse to the membrane. The model was then re- 
formulated into a Gibbs formulation, where the activity 
entered the force-balance equation of the membrane sur- 
face directly, instead of being part of the force-balance 
in the bulk fluids surrounding the membrane. An ex- 
pression for the active contribution to the force-balance 

in the Gibbs formulation, f^^^ = {r^ft^'lB + r„"a,ct»T') , 
was given through Equations (|33|l and (|40|l in the form of 
a moment expansion of the active force distribution, i^act- 
This expression provides a way to study the consequences 
of the force-dipole model of ^ for non-planar membrane 
shapes, opening up the possibility of testing the model in 
a wider class of experiments than theoretical predictions 
restricted to planar membranes would allow for. Some 
consequences of the activity, however, become immedi- 
ately apparent in the Gibbs formulation. One is that the 
dipole moment of the active force distribution tTdip gives 
rise to a modification of the tension of the membrane. 
The prominent effect of the next order in the moment 
expansion, the quadrupole moment Q, is that it modi- 
fies the bending moments of the membrane. The term 
in the force-balance of the membrane that this results 
in couples the dynamics of the membrane shape to the 
protein concentrations, and this term, if the force-dipole 
model is a proper explanation for the activity induced 
effect observed in the micropipette experiments 0, |3, 
induces a magnification of primarily the bending rigid- 
ity dominated part of fluctuation spectrum Q, i.e. the 
q"^ behavior of Equation (|53|) . It should be noted that 
this modiflcation of the fluctuation spectrum does not 
correspond simply to a renormalized bending rigidity or 
a modification of another coupling parameter from the 
equilibrium free energy. The reason behind this is that 
the activity does not produce symmetric couplings be- 
tween the dynamics of the membrane shape and the pro- 
tein densities. The coupling constants in the free energy 
produce such symmetric couplings because of Maxwell 
relations and Onsager's reciprocal relations. 

The results obtained in this paper are for a model of 



the activity that generalizes the force-dipole model of |j] . 
If one wants to reduce the equations describing the mem- 
brane dynamics here to those of 0| one first of all needs to 
specialize to the force distribution given by Equation 115|) 
|26l |. Secondly, one has to go to the planar case. Thirdly, 
thermal noise should be included in the Navier-Stokes 
equation. Equation (|12|) or H19|l depending on which for- 
mulation one begins with. Note that non-thermal noise 
due to the activity is not included in 0|. See 
for ways to include this. Fourthly, the bound- 
ary condition of Equation H18|) should be specified. In 
effects due to permeation of the membrane are in- 
cluded at the outset, but it is later argued that such 
effects are negligible for undulations with wavelengths 
in the regime relevant to experiments. Since bulk fiuid 
motion in the tangential direction of the membrane de- 
couples from the motion of the membrane shape in the 
planar case (see Equation ljA13l) here) it is enough to 
choose n ■ U± — n ■ dR/dt to obtain the equations of 
that govern the membrane motion in the experimen- 
tally relevant regime. Fifthly, diffusion equations for the 
proteins in the membrane should be set up. The hydro- 
dynamic attraction proposed in section IVII here can be 
included simply as a decreased diffusion constant of the 
active proteins. 

Because of the simplified bulk hydrodynamics in the 
Gibbs formulation, this formulation will constitute a 
much simpler starting point than the semi-microscopic 
formulation for further studies of long- wavelength conse- 
quences of the force-dipole model. This will be applied in 
another paper jlJl, where the fluctuation spectrum of a 
quasi-spherical vesicle will be derived for the force-dipole 
model. 
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APPENDIX A: SOLUTION OF THE 
SEMI-MICROSCOPIC FORMULATION IN THE 
PLANAR CASE 

In this appendix we will find the solution of the semi- 
microscopic formulation. Equation (|12l) . to first order in 
deviations from a planar shape situated at z = in a 
Cartesian coordinate system (x,y,z). 



1. Hydrodynamics without active forces 

We need the solution without the active contribution, 
such that we have a reference when we wish to find the 
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modification due to the activity. To find this solution we 
will Fourier transform in the plane of the membrane 

v±{x,y,z) - j ^e-'(«--+«»^)t;±,,(z) , (Al) 

and decompose the velocity into its z, longitudinal and 
transverse components as 



V±-q ^'"Lq~^ + V±,q^ + AJ 



(A2) 



where q = {qx,qy,Q) = qq and f is a unit vector perpen- 
dicular to both of the unit vectors z and q. The solution 
of H19|) can then be written 

= [(l + g|z|)4^,U.o + *9^4,,U=o] , (A3) 

[(l-g|z|)4,,U=o + *9^4,,U=o] , (A4) 

(A5) 

l^l,qU=o + iv'±J.=o ■ (A6) 



4,, 


= e^9l^ 












= 2r]qe' 



We will also find the forces of the bulk fluids on 
the membrane defined in Equation H24|l . To this end we 
need the stress tensors 



The important components at z = are 



(A7) 



2" •i"±.q 



e ""'gwj X^o , (A8) 



T2v 
TV 



(2^)2 
(Pq 



-i(qxX+qyy) t 



and they give the forces on the membrane 



Tt q 



2riqv^^ 



(A9) 
(AlO) 



(All) 
(A12) 



Kq ■ * = - mv±Jz=o , 

(T± • n)^ = T i27rf S (q,) S (qy) - 2r,qvl^U=o , 

(A13) 

where = ■ tat" is the projection of the forces onto 
the tangent space of the membrane. 



Expanding to first order in deviations e from the planar 
shape we get 

Fact,/,. ^F,,tdh {z(l + hV^e - e9,) - V||e} 5 (z - /i) , 

(A15) 

where V|| = f^da- Using this force in the Navier-Stokes 
Equation H12(l . together with the boundary condition 
that the velocity should vanish at z = 0, we get for 
\z\<\h\ 

pf =i^act./.,,e-«l''' [ T cosh(gz) - g/ie-^'^'] , (A16) 



"2^ 



^ ^act,Me-«l''l[-ze-«l^lq/i 



(g-i + |/i|)sinh(g|z|)-|z|cosh(g|z|)] , 



«±,9 = T ^Fact,/.,,e-'?l"l [{h ~ z) sinh(gz) + g/^ze-'l^l] 



(A17) 

-9|z|] 

(A18) 



and for \z\ > \h\ 
Pt =i^act,/,,,e-«l-l [ (sinh(g/.) - g/ie-^l^l) ] , (A19) 



1 ^ 



4.9 =^i^act,/.,ge-'?l^l [ - z(g/ie-''l"l - sinh((?/i)) 



act,/i,g^ 

± (g^"'^ sinh((7/i) — ft, cosh((j/i))] 



(A20) 



2v 



where 



4,, = T :^i^act,/,,,e"«l^l [(/i - z) sinh(g/i) + gfte-^l'^U] 

(A21) 

Fact,/!. = Fi^ctdh (l + W|e) , (A22) 

and it is assumed that h and z are positive for the 
part and negative for the "—"-part. For Vj. ^ we simply 
have Vj.^ = everywhere. 

We can write Equations (|A16p to HA2ip as 



4,g =e-'''' [(1 + 91^1) v'±,qU=o + iqzv'±J.=o] 

2, extra 



■ V. 



±,1 



(A23) 



•A,q =e [(1 - gf|z|) v^±Jz=o + ig2w±,«U=o] 



+ v'. 



extra 
±,9 ' 



p ^2r]qe' 



±, extra 
Pq 



(A24) 

:tra 

(A25) 



2. Hydrodynamics with active forces 

We will consider the contribution to Equation p2|) 
from a single "layer" of active force 



Fact,/, ^dh dA Factn(5' {r-iR + hn)) . (A14) 



where 

Vl,ql=0 



Fact, /!.,<; (g ^ sinh((7/i) - ft cosh(g/i)) , 

(A26) 

Ti^Fe.ct,h.qhsmh{qh) , (A27) 
2/7 
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and the extra part is zero for \z\ > \h\ and 

^act, h,q 



± 



2r] 



[h — z) cosh((7(z — K)) 



V 



extra 



+ q ^ sinh(g(z — h))\ 

Fact ,h.q 



P. 



±, extra 



2-q 



(A28) 

[{h- z)sinh{q{z-h))] , (A29) 



T F.^ct,h,q cosh{q{z - h)) , 



(A30) 



when l^l < \h\. This extra part gives a force on the 
membrane which is 

Tt^--''^ ■ z ^ F,,t,h,q{cos\v{qh) ~ qhsmh{qh)) , (A31) 



-,±,oxtia ^ ^ ^ -iFi,ct,h,qqhcosh{qh) . 



(A32) 



Equations l)A23(l to (|A25|I should be compared with 
Equations (|A3|) to (|A6|) . Doing this we see that we get 
the same behavior in the bulk for \z\ > \h\ if remove 
the active force from the Navier-Stokes equation and in- 
stead use the boundary conditions in Equations (|A26|I 
and HA27p . To second order in h these boundary condi- 
tions are 



'±,q\z 



=0 = 0, 



v±^q\z=o = Tiq-^F^ctji,qdhh^ 

Zl] 



(A33) 
(A34) 



which agrees with the boundary condition in the Gibbs 
formulation, Equation 144() . 



If the Gibbs formulation is correct, then Equations 
(|A31|I and (|A32(I should match the active force f^^^ from 
section ^3 Expanding Equations HA31|) and (jA32|) to 
second order in h we find 



+ , extra 
Q 



extra , /yi—, extra \ 
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) • z —Fact,h,qdh + 2Fi,ct,h,q=adh hHg 
- \q^F,,t^h.gdhh^ , (A35) 
iqFact,h,qdh h 



2iqHqFs,ctM,q=odhh'^ 



(A36) 



and again agreement with the Gibbs formulation. Equa- 
tions (|33|l and H4UI) . is found when it is recalled that the 
monopole moment of the force i^act vanishes. 



[1] B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, 

and P. Walter, Molecular Biology of the Cell (Garland, 

New York, 2002), 4th ed. 
[2] J.-B. Manneville, P. Bassereau, D. Levy, and J. Prost, 

Phys. Rev. Lett. 82, 4356 (1999). 
[3] P. Girard, J. Prost, and P. Bassereau, Phys. Rev. Lett. 

94, 088102 (2005). 
[4] J.-B. Manneville, P. Bassereau, S. Ramaswamy, and 

J. Prost, Phys. Rev. E 64, 021908 (2001). 
[5] N. Gov, Phys. Rev. Lett. 93, 268104 (2004). 
[6] J. Pecreaux, H.-G. Dobereiner, J. Prost, J.-F. Joanny, 

and P. Bassereau, Eur. Phys. J. E 13, 277 (2004). 
[7] U. Seifert, Adv. Phys. 46, 13 (1997). 
[8] M. A. Lomholt, P. L. Hansen, and L. Miao, Eur. Phys. 

J. E 16, 439 (2005). 
[9] J. W. Gibbs, The Scientific Papers of J. Willard Gibbs, 

Vol. 1: Thermodynamics (Dover, New York, 1961). 
[10] N. Kahya, D. A. Wiersma, B. Poolman, and D. Hoekstra, 

J. Biol. Chem. 277, 39304 (2002). 
[11] M. A. Lomholt, Phys. Rev. E 73, 061914 (2006). 
[12] S. Sankararaman, G. I. Menon, and P. B. S. Kumar, 

Phys. Rev. E 66, 031914 (2002). 
[13] W. Cai and T. C. Lubensky Phys. Rev. E 52, 4251 

(1995). 

[14] L. Miao, M. A. Lomholt, and J. Kleis, Eur. Phys. J. E 9, 
143 (2002). 

[15] M. A. Lomholt and L. Miao, J. Phys. A: Math. Gen. 39, 
10323 (2006). 

[16] Y. Kimura, D. G. Vassylyev, A. Miyazawa, A. Kidera, 



M. Matsushima, K. Mitsuoka, K. Murata, T. Hirai, and 

Y. Fujiyoshi, Nature 389, 206 (1997). 
[17] P. A. Kralchevsky, J. C. Eriksson, and S. Ljunggren, Adv. 

Coll. and Interface. Sci. 48, 19 (1994). 
[18] R. Capovilla and J. Guven, J. Phys. A: Math, and Gen. 

35, 6233 (2002). 
[19] J. R. Henriksen and J. H. Ipsen, Eur. Phys. J. E 14, 149 

(2004). 

[20] N. Dan, A. Berman, P. Pincus, and S. A. Safran, J. Phys. 

II 4, 1713 (1994). 
[21] T. Gulik-Krzywicki, M. Seigneuret, and J. L. Rigaud, J. 

Biol. Chem. 262, 15580 (1987). 
[22] P. Lenz, J.-F. Joanny, F. Jiilicher, and J. Prost, Eur. 

Phys. J. E 13, 379 (2004). 
[23] H.-G. Dobereiner and U. Seifert, Europhys. Lett. 36, 325 

(1996). 

[24] As suggested by a referee of this article, the experiment 
might be made even more sensitive to the protein activity 
if the dynamics of the bistability of the shape is studied 
at the point of the pro late- oblate transition in the way 
that it is studied in |23 | for the equilibrium situation. 

[25] For actual estimates of Reynolds numbers for membrane 
hydrodynamics the reader is referred to [3] 

[26] The moment expansion employed in this paper is not 
a restriction in comparison with Q, since the same ex- 
pansion is performed in Q (without it being mentioned 
explicitly though). 



